Cosmological Models with Shear and Rotation 
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ABSTRACT 

Cosmological models involving shear and rotation are considered, first in 
the General Relativistic and then in the Newtonian framework with the aim 
of investigating singularities in them by using numerical and analytical tech- 
niques. The dynamics of these rotating models are studied. It is shown that 
singularities are unavoidable in such models and that the centrifugal force 
arising due to rotation can never overcome the gravitational and shearing 
force over a length of time. 



I. INTRODUCTION 

Cosmological models with rotation have drawn attention in the past for 
various reasons. Godel first proposed a rotating universe model, with the 
aim of demonstrating that general relativity does not incorporate Mach's 
principle. Another reason for proposing the rotating models has been with 
the view to avoid a singularity in the universe. This approach was pioneered 
by Heckmann and Schucking. The basis for this line of thought has been 
that the centrifugal force arising due to rotation would prevent a collapse of 
the universe. Finally, since all astrophysical systems have rotation in them 
so it has been cojectured that there might be rotation in the universe too. 

The present work will deal with the investigation of singularities in a 
special class of metrics namely the Heckmann Schucking metrics in the Rela- 
tivistic and Newtonian framework. Although in general relativity a singularity 
is unavoidable, as shown by Hawking and Penrose, the study of the dynam- 
ics of rotating models does yield an idea of how far rotation may be able 
to prevent a collapse. The centrifugal force may prevent collapse along the 
axis but the collapse may nonetheless take place perpendicular to the axis 
of rotation. The first part of the work deals with the relativistic Heckmann 
Schucking model while in the second part the Newtonian analogues of rela- 
tivistic models will be looked at. 



II. THE RELATIVISTIC HECKMANN SCHt/CKING METRIC 

The most common approach to describe the cosmology of the universe is 
on the basis of the assumption of homogenity and isotropy of the universe 
which leads to the Robertson - Walker metric 

ds 2 = dt 2 - a 2 (t)[dr 2 /(l - kr 2 ) + r 2 (d6 2 + stn6 2 d(j) 2 )}, (1) 

where k=l,0,-l. 

This line element when coupled with the field equations, 

RI U, _ \g^ R = _ kT ^ _ A^, (2) 

can give the value of a(t) where T^ v is given by 

where p and u are the matter and radiation densities and p is the pressure. 

This form of the solution gives rise to a universe with a singular ori- 
gin. Can the singularity be averted if one considers a modification in the 
assumptions governing the choice of the line element? The postulate regard- 
ing isotropy of the universe is replaced by anisotropy which may be due to 
rotation and shear. The world lines of the galaxies are still geodesies along 
which the world time is measured but the t=constant surfaces are no longer 
orthogonal to the geodesies. This gives rise to cross product terms of the 
type g^dx^dt^fi = 1, 2, 3) in the metric. 

One such type of metric was proposed by Godel (1949), 

ds 2 = dt 2 + 2e Xl dtdx 2 - (dx 1 ) 2 + -e 2xi (dx 2 ) 2 - (dx 3 ) 2 . (4) 

The solution of the field equations with this metric indicate a rotating 
and stationary universe with the angular velocity of the various components 
given by 

^ = ~^-=av\K- (5) 

The non-vanishing of u 3 shows a rotation about the x 3 = constant axis. 

The failure of this model to account for the red shift of galaxies because of 
its stationary character prompted Heckmann and Schiicking (1958) to give a 



generalisation of the Godel's metric to obtain a universe with a non-singular 
origin and also to account for the red shift by making this model non-static. 
They hoped to find finite oscillating universes with maximum and minimum 
radii. 

One such model having Godel's model as a special case, has the line 
element 
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ds 2 = dt 2 +2e xl dtdx 2 -c ll (t)(dx 1 ) 2 -2c 12 (t)e xl dx 1 dx 2 +ac 11 (t)e 2xi (dx 2 ) 2 -S 2 (t)(dx 3 ) 

(6) 
which when solved with the field equations gives the following set of differ- 
ential equations : 
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Heckmann and Schiicking did not solve these equations which are too 
complicated for an analytical solution. Narlikar (1960) showed, however, 
that these models would not fulfill their original purpose. Later the interest 
in these models waned as the singularity theorems gained wide acceptance. 
Nevertheless we return to the Heckmann Schucking model to investigate 
its dynamical behaviour by numerical methods, to see how the singularity 
actually develops. 

The above equations were thus solved numerically with the initial values 
taken in correspondence with the Godel (1949) metric except for the intro- 
duction of a ciiterm to introduce a non static character in the metric. Thus, 



cii(0) = 1, 



(13) 



ci 2 (0)=0, (14) 

A = ~ (15) 

and different values of c'u and a were taken to highlight different cases. The 
parameters were varied as follows: 

ch = 0.1, 0.01, 0.001 (16) 

« = i-i (17) 

in correspondence with the two seperate cases pointed out by Narlikar (1960). 

The results for a = \ as shown in figure (1) shows a collapse in R and 
figure (2) reveals the small oscillations R makes before the collapse. Figure 
(3) shows a continuously shrinking x 3 dimension. The coefficient of (dx 3 ) 2 
in the metric goes to zero linearly showing a one dimesional singularity. R, 
which shows the extent of the universe in the other two dimensions also goes 
to zero[Figure (1)] pointing to a singularity in the volume after a period of 
oscillations since V oc R 2 S. The equations also imply a rapidly increasing 
angular velocity by consequence of uj — -h^ as pointed out by Godel (1950) 
in one of his discussions of the properties of the metric. 

The oscillations are sustained for a longer period of time if the initial rate 
of change of cnwith respect to time is large as this delays R from falling to 
zero as a consequence of 

R 2 = en — ttC^ — c^ 2 , (18) 

becoming smaller and smaller more slowly. 

It can also be seen from figures (4) and (5) that for the case a — — | 
oscillations are not possible in the Godel like universe and for the same 
other parameters the universe shrinks to a two dimensional singularity with 
a rapidly increasing SfFigure (5)], the coefficient of (dx 3 ) 2 and a decreasing 
i?[Figure (4)]. 

The volume of the universe given by V oc R 2 S also goes to zero suggesting 
that a singularity results in this case too. The results for a = — \ are in 
agreement with the analytical argument of Narlikar (1960) where he has 
shown that for a — — \ the universe does not oscillate between finite limits 
by considering (7), (8), and (9) coupled with the signature conditions. This 
form of a universe leads to a situation similar to the isotropic case. 



Another result which may be directly seen from the numerical simulation 
is that the universe for a = — \ is inherently unstable and collapses about 
an order of magnitude faster than for a = |, all other parameters being 
comparable. 



III. NEWTONIAN COSMOLOGIES WITH ROTATION 

One of the first works done on developing cosmology on Newtonian terms 
was by Milne and McCrea (1934). Their work was based on the concept of 
isotorpy and homogenity in the universe and they showed that models anal- 
ogous to the Friedmann Robertson Walker model could be obtained by using 
the Newtonian theory too. Heckmann and Schucking (1955) formulated the 
Newtonian equations of cosmology by solving the following three equations 
used to describe the general behaviour of a homogenous universe, but doing 
away with the assumption of isotropy: 

P- + div.w = 0, (19) 

v + (v.V)v = -0, (20) 



V 2 + A = AupG. (21) 

The Newtonian equations of cosmology formulated by them include both 
shear and rotation terms. Thus while rotation prevents the universe from 
collapsing, shear has the opposite effect. This is analogous to the general 
relativistic result derived by Raychaudhuri (1955). They tried to avoid sin- 
gularity by setting the shear equal to zero which is possible in the Newtonian 
framework. They were able to show that it is possible to avoid singularities 
because the rotation term dominates as R — > , hence preventing a collapse. 

This freedom to take the shear equal to zero in Newtonian models does not 
exist in General Relativistic models as was shown by Ellis (1967). Hence, it 
is not possible to avoid singularities in the General Relativistic models which 
give a more accurate description than their Newtonian analogs, in accordance 
with the Hawking-Penrose theorem (1975). 



However in 1963, Narlikar formulated the gravitational force approach 
which was based on the inverse square law rather than the Poisson equation, 

V 2 = AnGp. (22) 

This puts a greater restriction on <fi and implies that even though it might 
be possible to put the shear equal to zero initially, the time dependence of 
the shear has the effect of it acquiring some finite value in the course of time. 
Thus it is not obvious that singularity can be prevented. 

By the axioms of Newtonian mechanics we have a Euclidean space with 
rectangular co-ordinates x^/j, = 1, 2, 3) and an even-flowing uniform time, t. 
Thus at any time t the universe will present the same large scale view to all 
observers whose motion is idealized as the streaming of an ideal fluid. The 
velocity-distance relation can be written in the form, 

Vy. = H^Xv, (23) 

where H^ v is a function oft only, which appears in analogy with the Hubble's 
constant in the Robertson- Walker cosmology, but is direction dependent in 
accordance with the assumption of anisotropy. After writing the solution of 
(23) in the form, 

x» = a^{t)x° u , (24) 

A = de£||a M „|| and solving (19) and (20), Narlikar (1963) arrived at the 
following equation of motion, 

2 ~ 3 [ ' 

After replacing the time t by a dimesionless co-ordinate, 

r = (^)lt (26) 

and continously differentiating (25), Narlikar (1963) was able to reduce the 
six equations of motion to a single fourth order equation: 

A 2 A"" + 7 A A" - 4 A' 2 + 9 A = 0. (27) 

This fourth order nonlinear differential equation may be reduced by the fol- 
lowing transformations, 

A = F\ (^) 2 = X(F), F = e u ,^ = Y(X), (28) 



X = F' 2 ,Y = 2FF", (29) 

to the following second order equation, 

XY 2 (^-)+XY& 2 + (X + ^-)Y%+Y 2 -2XY+7Y-2X + 9 = 0. (30) 
dj\ dJi 2 dX 

The above equation was solved numerically by Narlikar (1963) for various 
initial conditions and all of them did show a singularity. Davidson and Evans 
(1973) investigated (25) further and all their numerical and analytical results 
also show a singularity 

But, since numerical analysis cannot be exhaustive of all initial conditions, 
(30) was analysed asymptotically to arrive at a more general result. We take 
the form of the asymptote as, 

Y = mX + K (31) 

Substituting Y in (30) and equating the coefficients of X 2 and X to zero 
we find the following asymptotes to the curves 

Y = K u (32) 

Y = -IX + K 2 , (33) 

Y = ^X - 3. (34) 

Here K\ and K 2 are constants and (34) is an exact solution and was also 
found by Narlikar (1963). The curves with asymptotes Y = —2X + K\ were 
also investigated by Narlikar (1963), and were found to go asymptotically as 

Y = —2X as is shown in figure (6). However the curves with asymptotes 

Y — K\ were missed in that work. 

Substituting the value of Y from (29) in (32), we get, 

FF" = k. (35) 

Solving the above equation by integrating leads to, 

dF 

C\dr. (36) 



VT^F 

Substituting x 2 = InF, we can get the equation in the form of the integral 
of the error function, 



-2 



e x dx = c 2 t, (37) 

where c 2 is some constant independent of A and r. 

It may be seen from (37) that r has a range of values from — oo to oo 
being the time co-ordinate. By readjusting the time scale we can always 
arrange that at some point in the range of r, A must become equal to zero 
which will imply a singularity 

The two other asymptotes of the curve have already been investigated 
numerically by Narlikar (1963) and Davidson and Evans (1973) and were 
shown to be singular. The asymptote Y— constant does indeed correspond 
to a few numerical solutions which were compiled by Davidson and Evans 
(1973). 

Another way of looking at the singular nature of the solution is by inves- 
tigating the case when asymptotically Y = —2X. Then we can take, 

Y = -IX + c= -2XZ(X) (38) 

where the function Z(X) — > 1 as X — > oo. Thus we get, 

Z{X) = 1 - ^ (39) 

where c is a constant independent of X. 

Substituting these values of Y and its first and second derivatives in (30) 
we see that the quadratic terms of (30) exactly cancel to zero. We are then 
left with 7Y — 2X + 9. For this to tend to zero asymptotically we arrive at 
the condition, 

- 16X + 7c + 9 = (40) 

Substituting the value of X from (29) and putting it in terms of A from 
(28) we have, 

|£ + 7c + 9 = (41) 

Thus, 

Ao -9 - 7c 

^ = ea*[(-^)(T 2 -7i)] (42) 

It is clearly seen from (42) that for c> — | A2 — *• for t 2 »t\ and for c< 
-f , A 2 -»• for r 2 « Ti. 



Since the time scale can be arbitrarily adjusted , we will get a singularity 
in one of the cases mentioned above. 

The conclusions arising from the fact that the volume element of the 
universe, A goes to zero in all possible cases and taken into consideration 
implies a singular solution to the generic Newtonian cosmological model with 
shear and rotation. The results show that in an expanding universe the 
outward force arising due to rotation can never be enough to overcome the 
combined force due to shear and gravitational pull. It may also be seen that 
the shear term cannot be taken equal to zero over the entire time scale if the 
expanding model of the universe is considered. 



IV. CONCLUSION 

Though singularities are unavoidable as shown by Hawking and Penrose 
(1975) in the general singularity theorem, nonetheless it is interesting to ob- 
serve in these models how we arrive at singularities in the course of time. 
This result holds equally well for the Newtonian models with shear and ro- 
tation as they have been shown to be singular in the history or future of the 
universe. The centrifugal force arising due to the rotation of the universe can 
never adequately combat the gravitational and shearing force. 
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Figure 1: Variation of R with time with a = 0.5. The figure shows the 
behaviour of R with time for a positive a and reveals a behaviour culminating 
in a collapse. 
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Figure 2: Variation of R with time with a = 0.5. The figure shows 
that during the initial time the behaviour of R for a positive a reveals an 
oscillating pattern. 
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Figure 3: Variation of S with time with a = 0.5. The figure shows the 
behaviour of the coefficient of the x 3 dimension, for a positive a, which is 
seen to go to zero linearly with time. 
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Figure 4: Variation of R with time with a = —0.5. The figure shows 
that for a negative a, R goes to zero monotonically without executing any 
oscillations about an order of magnitude faster than the case with a positive 
a. 
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Figure 5: Variation of S with time with a = —0.5. The figure shows 
that for a negative a, the coefficient of the x 3 dimension increases with time 
first linearly and then very rapidly. 
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Figure 6: Variation of X with Y for some different initial conditions. 

The curves show that the curves for Y > can never be joined to the curves 
for Y < and that Y eventually decreases with X corresponding to universes 
where rotataion has been ineffective in preventing a singular state. All curves 
in the figure asymptotically go as Y = —2X. 
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